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SECTIOr  I 
INTRODUCTION 


• It  of  c - min  chemical  species  into  the  upper  atmosphere  results 
' . iisplay  the  resonance  electronic-vibrational-rotal 

.jtrum  of  tn-  released  species.  Such  spectra  are  seen  in  rocket  releases 
■n.icals  f or  upper  atmospheric  studies  and  upon  reentry  into  the  atrnos- 
Llit<  . d s.  ( : particular  interest  in  this 

it  ion  i.  observed  spectra  of  certain  metallic  oxides  arid  air  diatomic 
i or i at  /: ic  r;  icies.  iron:  band  intensity  distribution  of  the  spectra,  and 
v.'ledre  of  the  f -values  for  electronic  and  vibrational  transitions,  the 
local  conditions  of  the  atmosphere  can  be  determined  (Ref.  1). 

I & 

NJ 

iTesont  theoretical  efforts  which  are  directed  toward  a more  complete 
■•-vi  realistic  analysis  of  the  transport  equations  governing  atmospheric 
relaxation  and  the  propagation  of  artificial  disturbances  require  detailed 
inf  rmation  of  thermal  opacities  and  LWIR  absorption  in  regions  of  temperature 
n.vi  pressure  where  molecular  effects  are  important  (Rofg.  ■ 'i± ■ and-— t) . Although 
.various  experimental  techniques  have  been  employed  for  both  atomic  and 
molecular  systems,  theoretical  studies  have  been  largely  confined  to  an  analy- 
sis of  the  properties  (bound-bound,  bound-free  and  free-free)  of  atomic  sys- 
tems 1 and  5-).  This  has  been  due  in  large  part  to  the  unavailability 

of  reliable  wavefunctions  for  diatomic  molecular  systems,  and  particularly 
for  excited  states  or  states  of  open-shell  structures.  More  recently,  fftgfs. 

~6*r9’)  reliable  theoretical  procedures  have  been  prescribed  for  such  systems 
which  have  resulted  in  the  development  of  practical  computational  programs.  __ 


The  theoretical  analysis  of  atmospheric  reactions  requires  the  knowledge 
of  the  electronic  structure  of  atoms,  ions  and  small  molecular  clusters  of 
nitrogen  and  oxygen  and  the  interaction  of  water  or  other  small  molecules 
with  these  clusters.  Knowledge  of  the  chemistry  of  metal  oxide  species,  which 
might  be  present  in  a contaminated  atmosphere,  is  also  desired.  Because  of 
the  computational  complexity  for  systems  with  large  numbers  of  electrons,  tra- 
ditional ab  initio  theoretical  methods  are  difficult  and  expensive  to  apply. 
However,  ab  initio  computational  programs  based  on  Slater-type  orbitals  and 
incorporating  both  SCF  and  configuration  interaction  analysis  have  been 
developed  and  are  available  for  studies  to  chemical  accuracy  for  systems  such 
as  Og,  N2,NO  and  their  corresponding  molecular  ions.  Semiempirical  methods, 
such  as  extended  Huckel  or  CNDO  theories  may  be  applied  to  systems  containing 
many  more  electrons  or  nuclear  centers  but  they  are  dependent  on  the  parameter 
ization  of  many-center  integrals,  which  is  difficult  and  often  ambiguous  for 
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systems  containing  many  different  atoms.  However,  there  is  a local  potential 
model  for  the  electronic  structure  of  atoms,  molecules,  clusters  and  solids, 
(Ref.  10),  which  is  relatively  simple  to  apply  to  systems  containing  heavy 
atoms  and  large  numbers  of  electrons.  This  is  a density  functional  method 
which  has  been  used  for  many  years  in  solid  state  theory  for  the  calculation 
of  the  energy  band  structure  of  solids  (Ref.  11),  and  has  been  shown  to  give 
accurate  results  for  heavy  atom  solid-state  systems,  both  for  the  one-electron 
spectrum  (Ref.  12)  and  the  cohesive  energy  (Ref.  13).  In  the  case  of  an 
infinite  perfect  crystal,  the  density  functional  method  is  usually  used  in 
conjunction  with  the  augmented  plane  wave  (APW)  or  the  Korringa-Kohn-Rostoker 
(KKR)  methods  of  solution  to  the  one-electron  Schrodinger  equations.  Until 
recently,  it  has  not  been  practical  to  apply  this  model  to  finite  systems 
other  than  atoms  (Ref.  14).  However,  in  the  last  several  years,  the  develop- 
ment of  "multiple-scattering"  techniques  (Refs.  15  and  l6)  for  the  solution 
of  the  equations  in  a multicenter  finite  system  has  allowed  the  density  func- 
tional model  to  be  applied  to  a large  number  of  molecules  (Refs.  17  through 
20)  and  clusters  of  atoms  simulating  localized  excitations  in  solids  (Refs. 

21  and  22)  which  cannot  be  described  by  an  energy  band  model.  In  its  usual 
form,  certain  assumptions  are  made  concerning  the  modecular  potential  in  order 
to  make  the  method  computationally  practical.  These  lead  to  the  so-called 
Y.a  computational  method.  During  the  present  research  program  under  Contract 
F19628-73-C-0300,  a thorough  study  of  the  effects  of  local  exchange  and  model 
potentials  was  carried  out  on  prototype  systems.  Substantial  improvements  in 
the  method  resulted  from  these  studies,  especially  in  the  reliability  of  the 
calculation  of  the  total  electronic  energy. 
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SECTION  II 

CURRENT  STATUS  OF  QUANTUM  MECHANICAL 
METHODS  FOR  DIATOMIC  SYSTEMS 


The  application  of  quantum  mechanical  methods  to  the  prediction  of 
electronic  structure  has  yielded  much  detailed  information  about  atomic  and 
molecular  properties  (Ref.  7).  Particularly  in  the  past  few  years,  the 
availability  of  high-speed  computers  with  large  storage  capacities  has  made 
it  possible  to  examine  both  atomic  and  molecular  systems  using  an  ab  initio 
approach,  wherein  no  empirical  parameters  are  employed  (Ref.  23).  Ab  initio 
calculations  for  diatomic  molecules  employ  a Hamiltonian  based  on  the  non- 
relativistic  electrostatic  interaction  of  the  nuclei  and  electrons,  and  a 
wavefunction  formed  by  anti symmetrizing  a suitable  many-electron  function  of 
spatial  and  spin  coordinates.  For  most  applications  it  is  also  necessary 
that  the  wavefunction  represent  a particular  spin  eigenstate  and  that  it  have 
appropriate  geometrical  symmetry.  Nearly  all  the  calculations  performed  to 
date  are  based  on  the  use  of  one-electron  orbitals  and  are  of  two  types: 
Hartree-Fock  or  configuration  interaction  (Ref.  8). 

Hartree-Fock  calculations  are  based  on  a single  assignment  of  electrons 
to  spatial  orbitals,  following  which  the  spatial  orbitals  are  optimized, 
usually  subject  to  certain  restrictions.  Almost  all  Hartree-Fock  calculations 
have  been  subject  to  the  assumption  that  the  diatomic  spatial  orbitals  are 
all  doubly  occupied,  as  nearly  as  possible,  and  are  all  of  definite  geomet- 
rical symmetry.  These  restrictions  define  the  conventional,  or  restricted, 
Hartree-Fock  (RHF)  method  (Refs.  24  and  25).  RHF  calculations  can  be  made 
with  relatively  large  Slater-type  orbital  (STO)  basis  sets  for  diatomic 
molecules  with  first  or  second -row  atoms,  and  the  results  are  convergent  in 
the  sense  that  they  sire  insensitive  to  basis  enlargement.  The  RHF  model  is 
adequate  to  give  a qualitatively  correct  description  of  the  electron  inter- 
action in  many  systems,  and  in  favorable  cases  can  yield  equilibrium  inter- 
atomic separations  and  force  constants.  However,  the  double-occupancy 
restriction  makes  the  RHF  method  inappropriate  in  a number  of  circumstances 
of  practical  interest.  In  particular,  it  cannot  provide  potential  curves 
for  molecules  dissociating  into  odd-electron  atoms  (e.g.,  NO  at  large  inter- 
nuclear  separation),  or  into  atoms  having  less  electron  pairing  than  the 
original  molecule  [e.g.,  On  ^g~  "*  0(^P)];  it  cannot  handle  excited  states 
having  unpaired  electrons  (e.g.,  the  states  of  O2  responsible  for  the 
Schumann-Runge  bands);  and,  in  general,  it  gives  misleading  results  for 
molecules  in  which  the  extent  of  electron  correlation  changes  with 
intemuclear  separation. 
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Confi . urat ion-interaction  (CIj  methods  have  the  c a:  ability  of  avoiding 
the  limi  tat  ions  f th<  RHP  calculations.  ! f conf  i g irai  Lorn  noi  r<  tricted 
doubly-occupii  l orbitals  are  incJ  . : 1,  a :,'l  can,  in  principle,  converge  ,c 
an  exact  wave!  : :tion  for  the  customary  Hamiltonian.  However,  mar:y  Cl  calcula- 
tions have  in  fact  been  based  on  a restriction  to  doubly-occupied  orbitals 
and  therefore  retain  many  of  the  disadvantages  of  the  RHF  method  (Ref.  8), 

The  use  of  general  Cl  formulations  involves  three  considerations,  all  of 
which  nave  been  satisfactorily  investigated:  the  choice  of  basis  orbitals, 

the  choice  of  configurations  (sets  of  orbital  assignments),  and  the  specific 
calculations  needed  to  make  vavefunctions  describing  pure  spin  states 

if,  6).  The  first  consideration  is  the  art  associated  with  quantum  mechan- 
ical electronic  struct’ire  calculations.  Many  methods  (iterative  MSO,  per- 
turbation selection,  first  order  Cl,  etc.)  have  been  advocated  for  the  optimum, 
choice  of  configurations.  T'nerri  are  no  firm  rules  at  present  and  the  optimum 
choice  is  a strong  function  of  the  insight  of  the  particular  research  investi- 
gator. The  last  consideration,  proper  s:  in  and  symmetry  projection,  has 
proved  difficult  to  implement,  but  computer  programs  have  been  developed  for 
linear  projection  algebra  at  this  Center,  and  the  Cl  method  has  been  found  of 
demonstrable  value  in  handling  excited  states  and  dissociation  processes 
which  cannot  be  treated  with  RHF  techniques. 

Either  of  the  above  described  methods  for  ab_  initio  calculations  reduces 
in  practice  to  a series  of  steps,  the  most  important  of  which  axe  the  evalua- 
tion of  molecular  integrals,  the  construction  of  matrix  elements  of  the  Hamil- 
tonian, and  the  optimization  of  molecular  orbitals  (RHF)  or  configuration 
coefficients  (Cl).  For  diatomic  molecules,  these  steps  are  all  comparable  in 
their  computing  time,  so  that  a point  has  been  reached  where  there  is  no 
longer  any  one  bottleneck  determining  computation  speed.  In  short,  the 
integral  evaluation  involves  the  use  of  ellipsoidal  coordinates  and  the  intro- 
duction of  the  Neumann  expansion  for  the  interelectronic  repulsion  potential 
(Ref.  27);  the  matrix  element  construction  depends  upon  an  analysis  of  the 
algebra  of  spin  eigenfunctions  (Ref.  28);  and  the  orbital  or  configuration 
optimization  can  be  carried  out  by  eigenvalue  techniques  (Refs.  29,  30).  All 
the  steps  have  by  now  become  relatively  standard,  and  can  be  performed 
efficiently  on  a computer  having  32,000  to  65,000  words  of  core  storage,  a 
cycle  time  in  the  microsecond  range,  and  several  hundred  thousand  words  of 
peripheral  storage. 

Both  the  RHF  and  Cl  methods  yield  electronic  wavefunctions  and  energies 
as  a function  of  the  internuclear  separation,  the  RHE  method  for  one  state, 
and  the  Cl  method  for  all  states  considered.  The  electronic  energies  can  be 
regarded  as  potential  curves,  from  which  may  be  deduced  equilibrium  inter- 
nuclear separations,  dissociation  energies,  and  constants  describing  vibra- 
tional and  rotational  motion  (including  anharmoric  and  rotation-vibration 
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effects ).  It  is  also  possible  to  solve  the  Schrodinger  equation  for  the 
motion  of  the  nuclei  subject  to  the  potential  curves,  to  obtain  vibrational 
wavefunctions  for  use  in  transition  probability  calculations.  The  electronic 
wave  functions  themselves  can  be  used  to  estimate  dipole  moments  of  individual 
electronic  states,  transition  moments  between  different  electronic  states,  and 
other  properties.  While  all  of  the  calculations  described  in  this  paragraph 
have  been  carried  out  on  some  systems,  the  unavailability  of  good  electronic 
wavefunctions  and  potential  curves  has  limited  actual  studies  of  most  of  these 
properties  to  a very  small  number  of  molecules. 

It  is  not  always  possible  or  practical  to  perform  an  ab  initio  or  mathe- 
matically rigorous  calculation  of  the  electronic  structure  of  a diatomic 
molecule  containing  a large  number  of  electrons  and  it  is  sometimes  useful 
to  resort  to  approximations  based  on  certain  physical  insights.  One  such 
approximation  is  the  density  functional  method  (Ref.  31)  as  applied  to 
molecules.  In  one  implementation  of  this  approach,  namely  the  model, 
the  one-electron  Schrodinger  equation  is  set  up  within  an  approximation  to  the 
true  potential  which  is  spherically  symmetric  within  spheres  surrounding  the 
several  nuclei,  constant  in  the  region  between  adjacent  spheres  and  spheri- 
cally symmetric  outside  an  outer  sphere  which  encloses  the  entire  collection 
of  nuclei.  The  resulting  equations  are  solved  by  a multiple-scattering  method, 
equivalent  to  the  KKR  method  (Ref.  32)  often  used  for  crystals.  From  the 
resulting  charge  densities,  onecan  compute  a total  energy,  using  a statistical 
approximation  for  the  exchange  correlations.  This  approximation  is  defined 
apart  from  a multiplicative  factor,  o,  hence  the  name  X^  method. 

The  method  has  been  very  successful  in  predicting  excitation  energies 
in  molecules  and  molecular  ions  but  has  met  with  only  limited  success  in  cal- 
c’ilation  of  total  energies.  The  basic  deficiency  appeared  to  lie  in  the 
muffin-tin  approximation  to  the  true  potential  rather  than  in  the  statistical 
treatment  of  the  correlation  energy.  However,  recent  studies  at  this  Center 
have  shown  that  when  certain  variational  constraints  are  applied  to  the 
wavefunctions,  significant  improvements  in  the  calculated  total  energies  can 
be  obtained.  At  this  time  the  method  appears  to  have  utility  for  constructing 
regions  of  potential  energy  surfaces. 
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SECTION  III 
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METHOD  OF  APPROACH 


A.  Quantum  Mechanical  Calculations 

Central  to  the  proposed  studies  are  the  actual  quantum-mechanical 
calculations  which  must  he  carried  o it  for  the  molecular  structures.  For 
added  clarity,  various  aspects  of  these  calculations  are  discussed  in 

individual  subsections . 


1.  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  near 
adequacy  of  a minimum  Slater- type-orbital  (STO)  basis  for  constructing 
molecular  wavefunctions  (Refs.  33  and  3M • This  means  inner-shell  and 
valence-shell  STO's  of  quantum  numbers  appropriate  to  the  atoms  (is,  2s,  2p 
for  C,  N,  0;  Is  for  H) . The  main  deficiency  of  the  minimum  basis  set  is  in 
its  inability  to  describe  polarization  of  n orbitals  in  atoms  adjacent  to  H 
atoms,  and  successful  calculations  usually  result  if  one  (or  a set)  of  p 
orbitals  is  supplied  for  each  H atom.  Values  of  the  screening  parameters  £ 
for  each  orbital  can  either  be  set  from  atomic  studies  or  optimized  in  the 
molecule;  the  latter  approach  is  indicated  for  studies  of  maximum  precision. 
When  high  chemical  accuracy  is  required,  as  for  the  detailed  studies  of  the 
ground  state  of  a system,  a more  extended  basis  should  be  used.  Double- 
zeta  plus  polarization  functions  or  optimized  MO's  usually  are  required. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavefunctions  with  all 
kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  states  can  adjust  to  an  optimum  form  at  each  geometrical  confor- 
mation and  for  each  state.  Except  when  well-defined  pairings  exist  for  as 
many  electrons  as  possible,  a single -configuration  study  (even  of  Hartree- 
Fock  quality)  will  be  inadequate. 


2.  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 
possess  definite  eigenvalues  of  the  spin  operator  S^  as  well  as  an  appro- 
priate geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by 
the  assignment  of  orbitals  to  each  configuration,  but  the  spin  state  must 
be  obtained  by  a constructive  or  projective  technique.  Formulas  have  been 
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developed  (Ref.  25)  for  projective  construction  of  spin  states  from 
orthogonal  orbitals,  and  programs  implementing  these  formulas  have  for 
several  years  been  in  routine  use  at  UTRC. 

One  of  the  least  widely  appreciated  aspects  of  the  spin-projection 
problem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes 
be  coupled  to  give  more  than  one  overall  state  of  given  S quantum  number. 

It  is  necessary  to  include  in  calculations  all  such  spin  couplings,  as  the 
optimum  coupling  will  continuously  change  with  changes  in  the  molecular 
conformation.  This  is  especially  important  in  describing  degenerate  or  near- 
degenerate excited  electronic  states. 


3.  Method  of  Ab  Initio  Calculation 


A spin-free  nonrelativistic  Hamiltonian  will  be  employed  in  the 
Born-Oppenheimer  approximation.  For  a molecule  containing  n electrons  and 
m nuclei,  the  approximation  leads  to  an  electronic  Hamiltonian  depending 
parametrically  on  the  nuclear  positions : 


J (R) 


Z a Z e 
+ R 


+ 


I 


i >i>  i 


(1) 


where  Zy.  and  % are  the  charge  and  position  of  nucleus  r^  is  the  position 
of  electron  i,  and  -??  is  the  Laplacian  operation  for  electron  i.  The 
Hamiltonian  H is  in  atomic  units  (energy  in  Hartrees , length  in  Bohrs). 

Electronic  wavefunctions  ty,  depending  parametrically  on  the  nuclear 
positions  are  made  to  be  optimum  approximations  to  solutions  of  the 
Schrodinger  equation 

^(RH'(R)=E(R)'HR)  (2) 

by  invoking  the  variational  principle 

8w(R,  = si£l5M i52SJ22i  (3i 

fV(R)V{R)dT 

The  integrations  in  Eq.  (3)  are  over  all  electronic  coordinates  and  the 
stationary  values  of  W are  approximations  to  the  energies  of  states  des- 
cribed by  the  corresponding  j[.  States  of  different  symmetries  are  studied 
by  restricting  to  the  appropriate  form,  and  excited  states  of  any  symmetry 
can  be  handled  by  simultaneously  determining  the  for  the  excited  state  and 
all  lower-lying  states  of  the  same  symmetry. 


The  specific  form  for  may  be  written 


'Hr)  = Z cM  (r)  (4) 

h 

where  each  is  referred  to  as  a configuration,  and  has  the  general 
structure 

.n  (iilR)oM  (5) 

where  each  • is  a spatial  orbital,  A is  the  antisymmetrizing  operator, 

0 ■ is  the  spin-projection  operator  for  spin  quantum  number  S,  and  cv,  is 
a product  of  a and  q one-electron  spin  functions  of  magnetic  quantum  number 
M.  TIo  requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbi- 
tal, no  Eqs.  (4)  and  (5)  can  describe  a completely  general  wavefunction. 

In  Hartree-Fock  calculations  is  restricted  to  a single  which  is 
assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 
orbitals  A . are  then  selected  to  be  the  linear  combinations  of  basis  orbi- 
tals  best  satisfying  Eq.  (3).  Writing 

V (6) 

The  a are  determined  by  solving  the  matrix  Hartree-Fock  equations 

l Sx„o„i  (each  X)  (7) 

where  e;  is  the  orbital  energy  of  . 
i 

The  Fock  operator  F^v  has  been  throughly  discussed  in  the  literature 
(Ref.  36)  and  depends  upon  one-  and  two-electron  molecular  integrals  and  upon 
the  asi.  This  makes  Eq.  (7)  nonlinear  and  it  is  therefore  solved  iteratively. 
UTRC  has  developed  programs  for  solving  Eq.  (7)  for  both  closed  and  open- 
shell  systems,  using  basis  sets  consisting  of  Slater- type  atomic  orbitals. 
Examples  of  their  use  are  in  the  literature  (Ref.  37). 

In  configuration  interaction  calculations,  the  summation  in  Eq.  (4) 
has  more  than  one  term,  and  the  c^  are  determined  by  imposing  Eq.  (3),  to 
obtain  the  secular  equation 

Z ws^„)  = 0 (each/i)  (8) 
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where  HM*  (R) ^ (R)'MR)^r 

S = j Vp  (R) 'i' t/  (R)  dT 

.quation  (7)  is  solved  by  matrix  diagonal! zation  using  either  a modified 
ens  method  (Ref.  2>)  or  a method  due  to  Shavitt  (Ref.  30 ) . 

Phi  atrix  ■ Lements  H and  S„  may  be  reduced  by  appropriate  erator 

u 

a’:  -jsra  to  the  forms 

€P  <(eM  0S  ? (r, , r)  J (r)  p n n,  (i,.r)S  (10) 

fP  (i,.R)  P n^jlrj.R^  (n) 

where  P is  a permutation  and  er  its  parity.  The  sum  is  over  all  permutations . 
0 5 P is  a "Sanibel  coefficient"  and  the  remaining  factors  are 

patia]  inte  ra  whicl  can  be  factored  into  one-  and  two-electron  integrals. 
Ef  the  ’2  . are  orthonormal,  Eqs.  (10)  and  (11)  become  more  tractable  and  thi 
H and  S.,,.  may  be  evaluated  by  explicit  methods  given  in  the  literature 
(Ref.  35) • iter  programs  have  been  developed  for  carrying  out  this  pro- 

cedure, and  they  have  been  used  for  problems  containing  up  to  bo  total  elec- 
trons, 10  unpaired  electrons,  and  several  thousand  configurations. 

The  01  studies  described  above  can  be  carried  out  for  any  orthonormal 
set  of  * ; for  which  the  molecular  Integrals  can  be  calculated.  Programs 
developed  by  UTP>C  made  specific  provision  for  the  choice  of  the  dpi  as 
Slater-type  atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree- 
fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 


h . Molecular  Integrals 

The  one-  and  two-electron  integrals  needed  for  the  above  described 
method  of  calculation  are  evaluated  for  STO's  by  methods  developed  by  the 
present  investigators  (Ref.  39)-  All  needed  computer  programs  have  been 
developed  and  fully  tested  at  UTRC. 


5.  Configuration  Selection 

Using  a minimum  basis  plus  polarization  set  of  one-electron  functions,  a 
typical  system  will  have  of  the  order  of  1(A  configurations  in  a full  Cl 
(that  resulting  from  all  possible  orbital  occupancies).  It  is  therefore 
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essential  to  identify  and  use  the  configurations  describing  the  significant 
part  of  the  wavefunction.  There  are  several  ways  to  accomplish  this  objective. 
First,  one  may  screen  atomic-orbital  occupa-cies  to  eliminate  those  with 
excessive  formal  charge.  Alternatively,  in  a molecular-orbital  framework 
one  may  eliminate  configurations  with  excessive  number  of  antibonding  orbi- 
tals. A third  possibility  is  to  carry  out  an  initial  screening  of  configura- 
tions, rejecting  those  whose  diagonal  energies  and  interaction  matrix  elements 
do  not  satisfy  significance  criteria.  Programs  to  sort  configurations  on  all 
the  above  criteria  are  available  at  UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (Ref.  40)  or  multiconfiguration 
SCF  (Ref.  4l)  concepts.  To  implement  the  natural-orbital  approach,  an  initi- 
al limited-CI  wave-function  is  transformed  to  natural-orbital  form,  and  the 
resulting  natural  orbitals  are  used  to  form  a new  Cl.  The  hoped-for  result 
is  a concentration  of  the  bulk  of  the  Cl  wavefunction  into  a smaller  number 
of  significant  terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome, 
but  in  principle  more  effective.  It  yields  the  optimum  orbital  choice  for 
a pre-selected  set  of  configurations.  This  approach  works  well  when  a small 
number  of  dominant  configurations  can  be  readily  identified. 

It  should  be  emphasized  that  the  problem  of  configuration  choice  is  not 
trivial,  and  represents  an  area  of  detailed  study  in  the  proposed  research. 

The  existence  of  this  problem  causes  integral  evaluation  to  be  far  from  a 
unique  limiting  factor  in  the  work. 


6.  Density  Functional  Approach  - X<*  Model 


The  Xqi  model  (Ref.  10)  for  the  electronic  structure  of  atoms,  molecules, 
clusters  and  solids  is  a local  potential  model  obtained  by  making  a simple 
approximation  to  the  exchange-correlation  energy.  If  we  assume  a nonrelati- 
vistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be  shown  that 
the  total  energy  E of  a system  can  be  written  exactly  (Ref.  42)  (in  atomic 


units)  as 
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This  expression  is  exact,  provided  the  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  Eq.  (12)  represents  the  kinetic 
and  electron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion 


energy.  The  sums  (u,,v)  are  over  all  the  nuclear  charges  in  the  system.  The 
third  term  is  the  electron-electron  repulsion  term,  which  represents  the 
classical  electrostatic  energy  of  the  charge  density  p interacting  with 
itself,  where 

pU)  =Z  n.u*(D  Uj(l)  (13) 


The  last  term  Exc  represents  the  exchange  correlation  energy  and  can  he 
expressed  formally  as 


E,c  z-rfp{,)  <*r,  i fief'll  d?r  . 


12 


(14) 


where  oxc  (l,2)  represents  the  exchange-correlation  hole  around  an  electron 
at  position  1.  In  the  exact  expression,  pxc  is  dependent  on  the  second-order 
density  matrix.  In  the  Hartree-Fock  approximation  Exc  is  the  exchange  energy, 
oxc  represents  the  Fermi  hole  due  to  the  exclusion  principle  and  depends  only 
on  the  first-order  density  matrix.  In  the  X&  method,  we  make  a simpler 
assumption  about  Pxc.  If  we  assume  that  the  exchange-correlation  hole  is 
centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown  that 
the  exchange-correlation  potential 

A / i r%  \ — * 

dri 


-I 


*C  J '12 

is  inversely  proportional  to  the  range  of  the  hole,  rg,  where  r„ 


(15) 

is  defined 
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Therefore,  in  the  X#  model,  the  potential  Uxc  is  proportional  to 
We  define  a scaling  parameter  a such  that 


p1/3M. 


u„a,l)=  - 
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( "bp  (J)/8tt) 


1/3 


(17) 


The  expression  in  Eq.  17  is  defined  so  that  a = 2/3  for  the  case  of  a free 
electron  gas  in  the  Hartree-Fock  model  (Ref.  43)  and  a = 1 for  the  potential 
originally  suggested  by  Slater  (Ref.  44) . A convenient  way  to  choose  this 
parameter  for  molecular  and  solid  state  applications  is  to  optimize  the 
solutions  to  the  Xa  equations  in  the  atomic  limit.  Schwarz  (Ref.  45)  has 
done  this  for  atoms  from  Z = 1 to  Z = 4l  and  found  values  between  2/3  and  1. 


In  the  "spin  polarized"  version  of  the  Xq-  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  only  with 
a potential  determined  by  the  charge  density  of  the  same  spin.  In  this  case 
the  contribution  to  the  total  energy  is  summed  over  the  two  spins,  s = i j, 


*c 


-r  i JVn  u.a,»(l)dfI 


(18) 
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where  the  potential  is  spin-dependent 

(l)=-  2 


U 


»a,, 


— ( 5p{  I )/4rr)'/3 


(19) 


and  o,.  is  the  charge  density  corresponding  to  electrons  of  spin  s.  The 
spin-polarized  Xa  niodel  is  useful  for  describing  atcrcs  and  molecules  with 
open-shell  configurations  and  crystals  which  are  ferromagnetic  or  anti- 
ferromagnetic. 


Once  one  has  made  the  approximation  to  the  total  energy  functional 
5 in  Eq.  (12),  then  the  rest  of  the  theory  follows  from  the  application  of  the 
variational  orinciole.  The  orbitals  u. 


».j  are  determined  by  demanding  that  E 
be  stationary  with  respect  to  variations  in  u. . This  leads  to  the  set  of  one- 
electron  'Xa  equations 
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(20) 


where  e-  is  the  one-electron  eigenvalue  associated  with  u^.  Since  0 (r)  is 


defined  in  terms  of  the  orbitals  u,-  , the  Eq.  (20)  must  be  solved  iteratively, 
•until  self-consistency  is  achieved.  Empirically,  if  one  takes  as  an  initial 
guess  that  o is  approximately  a sum  of  superimposed  atomic  charge  densities, 
then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor  of  2/3 


multiplying  the  potential  is  a result  of  the  nonlinear  dependence  of  Ev„  on 


o . This  also  has  as  a consequence  that  the  Xa  eigenvalues  e-  do  not  satisfy 
Koopmans ' theorem,  i.e.,  they  cannot  be  interpreted  as  ionization  energies. 
However  it  can  be  shown  that  the  e_.  are  partial  derivatives  of  the  total 
expression  of  Eq.  (12)  with  respect  to  the  occupation  number, 


€ . 
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(21) 


If  E were  a linear  function  of  n*,  then  Koopmans'  theorem  would  hold. 

However,  because  of  the  dominant  Coulomb  term,  E is  better  approximated  by 
a quardratic  function  in  n. . This  leads  to  the  "transition  state"  approxima- 
tion which  allows  one  to  equate  the  difference  in  total  energy  between  the 
state  (n^,  nj)  and  (n^  - 1,  + 1)  to  the  difference  in  the  one-electron 

energies  e j - e j_  calculated  in  the  state  (n^  - ■§■»  + |-).  The  error  in 

this  approximation  is  proportional  to  third-order  derivatives  of  E with 
respect  to  n^  and  n.,  which  are  usually  small  (Ref.  46).  The  main  advantage 
of  using  the  transition  state  rather  than  directly  comparing  the  total  energy 
values  is  computational  convenience,  especially  if  the  total  energies  are 
large  numbers  and  the  difference  is  small. 
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The  relationship  of  Ea.  (21) 
level"  for  the  ground  state.  This 


ni  under  the  condition  that  the 


also  implies  the  existence  of  a "Fermi 
can  be  seen  by  varying  E with  respect  to 
is  a constant,  i.e., 


sun  r n. 
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5 E-Xl  njJ  = ' 


(22) 


implies  *E/an^  = >,  where  \ is  a Lagrangian  multiplier.  This  implies  that 
the  total  energy  is  stationary  when  all  the  one-electron  energies  are  equal. 
However,  the  occupation  numbers  are  also  subject  to  the  restriction 


0 s n- 


This  leads  to  the  following  conditions  on  the  ground  state  occu- 


ation  numbers ; 


« <X  ) n.  = l 

> i 

«i  > X ) rij  =0 
c,  = X •)•  0<  n.  < i 


(23) 


In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with 
\ representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n^  are  either  0 or  1,  the  X® model  predicts , 
in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level.  In  particu- 
lar, this  will  occur  in  a system  (such  as  transition  metal  or  actinide  atom) 
which  has  more  than  one  open  shell. 


The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree-Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  Eq.  (20)  is  orbital-independent  and  local,  i.e.,  it  is  the  same 
for  all  electrons  (except  in  the  spin-polarized  Xa  theory)  and  is  a multipli- 
cative operator.  On  the  other  hand,  the  Hartree-Fock  potential  is  nonlocal, 
or  equivalently,  there  is  a different  local  potential  for  each  orbital. 

This  involves  a great  deal  more  computational  effort,  especially  for  systems 
described  by  a large  number  of  orbitals.  It  has  been  shown  (Ref.  47)  that 
the  Xa  orbitals  for  the  first  and  second  row  atoms  are  about  as  accurate  as 
a double-zeta  basis  set,  and  are  probably  better  for  larger  atoms  which 
involve  electrons  with  | s 2. 


2.  The  orbital-independent  Xa  potential  leads  to  a better  one-electron 
description  of  electronic  excitations  of  a system.  Both  the  unoccupied  (n^=0) 
and  occupied  (n^  = 1)  eigenfunctions  are  under  the  influence  of  the  same 
potential  resulting  from  the  other  N-l  electrons.  The  Hartree-Fock  virtual 
orbitals  see  a potential  characteristic  of  the  N occupied  orbitals,  and  there- 
fore are  not  as  suitable  for  describing  the  excited  states.  Actually, 
although  the  ground  state  virtual  eigne values  are  usually  a good  description 
of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
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potential  where  one-half  an  electron  has  been  removed  from  the  system  gives 
a much  better  first-order  picture  of  these  levels  (Ref.  48). 


3.  As  has  been  shown  by  Slater  (Ref.  49),  the  Xa  model  rigorously 
satisfies  both  the  virial  and  Heilman- Feynman  theorems,  independent  of  the 
value  of  the  parameter  a.  This  is  convenient  for  calculating  the  force  on 
a nucleus  directly  in  terms  of  a three-dimensional  integral,  rather  than  the 
six- dimensional  integrals  in  the  expression  for  the  total  energy  of  Eq.  (12). 


7.  Computational  Aspects  of  the  Xa  Method 

In  application  of  the  Xa  model  to  finite  molecular  systems,  there  are  two 
practical  aspects  of  the  calculations  which  must  be  considered.  The  first 
concerns  the  choice  of  the  integration  framework  for  describing  the  molecular 
wavefunctions  and  the  second  deals  with  the  choice  of  the  exchange  parameter, 
nt,  in  different  regions  of  space. 

In  computations  with  heteronuclear  molecules,  there  are  two  free  param- 
eters that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic  spheres 

of  integration  at  a given  internuclear  separation  and  the  value  of  the  exchange 
parameter  in  the  intersphere  region.  The  values  of  the  exchange  parameters 
are  usually  chosen  to  be  the  free  atom  values  inside  the  spheres. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a heteronuclear  diatomic  molecule  introduces  changes  in  the  total 
energy  that  can  be  large  on  a chemical  scale,  (~  1 eV) . A choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a good  esti- 
mate for  these  calculations.  Studies  made  at  UTRC  have  shown  that  at  any  given 
separation  the  total  energy  calculated  from  the  Xa  model  is  a minimum  at  the 
radii  ratio  where  the  spherically  averaged  potentials  from  the  two  atomic 
centers  are  equal  at  the  sphere  radius, 


This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic  species, 
the  long  range  tail  of  the  potential  must  go  like  -*2/R  from  one  ion  and  -2/R 
(in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separations,  the 
tails  of  the  potential  cannot  match  well.  However,  at  reasonable  separations, 
the  l/R  character  of  the  potential  does  not  invalidate  the  potential  match 
criterion  for  radii  selection.  This  match  for  the  atomic  potentials  is  applied 
to  the  self-consistent  potentials. 

The  second  free  parameter,  the  intersphere  exchange  coefficient,  is 
chosen  to  be  a weighted  average  of  the  atomic  exchange  parameters  from  the 
two  constituents.  At  small  internuclear  separations,  the  optimum  radius  for 
an  atomic  sphere  frequently  places  significant  amounts  of  charge  outside  that 
atomic  sphere  - charge  that  is  still  strongly  associated  with  its  original 
center  rather  than  being  transferred  to  the  other  center  or  associated  with 
the  molecular  binding  region.  To  best  account  for  these  cases  the  weighting 
coefficients  are  chosen  to  reflect  the  origin  of  the  charge  in  the  intersphere 
(or  outersphere  region), 
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where  (Qs^  - Qj°)  is  the  charge  lost  from  sphere  i relative  to  its  atomic 
value  (or  ionic  value)  Qj°  and  is  the  atomic  exchange  parameter  for  sphere 
i.  This  value  for  ^intersphere  calculated  dynamically  - it  is  updated 
after  each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  ex’s  for  small  atoms  vary  rapidly  with  z (and  with  ionicity) . The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
calculated  for  the  molecule  but  also  in  some  cases  affects  the  distribution  of 
charge  between  the  atomic  spheres  and  the  intersphere  region. 


B.  Transition  Probabilities 

The  electronic  and  vibrational -rotational  wavefunctions  of  a pair  of 
states  can  be  used  to  calculate  transition  probabilities.  If  two  molecular 
states  are  separated  in  energy  by  an  amount  AE^  = hcv  (h  = Planck's  constant, 
c = velocity  of  light,  v = frequency  in  wave  numbers),  the  semiclassical  theory 
of  radiation  (Refs.  50  and  51 ) yields  for  the  probability  of  a spontaneous 
transition  from  an  upper  state  n to  a lower  state  m 
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Here  is  the  Einstein  coefficient  for  spontaneous  transition  from  le  /el 
n -*  m,  gn  is  the  total  degeneracy  factor  for  the  upper  state 


9n  = (2-  Scr,A,)(2S'+  l)(2J'  + l) 


(25) 


and  is  the  total  strength  of  a component  line  is  a specific  state  of  polar- 
ization and  propagated  in  a fixed  direction.  A related  quantity  is  the  mean 
radiative  lifetime  of  state  n defined  by 
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-.he  summation  being  over  all  lower  levels  which  offer  allowed  connections, 
he  intensity  of  the  emitted  radiation  is 
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(27) 


where  K,  is  tne  number  density  in  the  upper  state  n.  This  analysis  assumes 
-hat  ail  degenerate  states  at  the  same  level  n are  equally  populated,  which 
will  be  true  for  isotropic  excitation.  The  total  line  strength  Snm  can  be 
written  as  the  square  of  the  transition  moment  summed  over  all  degenerate 
components  of  the  molecular  states  n and  m: 
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where  j and  i refer  to  all  quantum  numbers  associated  collectively  with  upper 
and  lower  electronic  states,  respectively. 

In  the  Born-Oppenheimer  approximation,  assuming  the  separability  of 
electronic  and  nuclear  motion,  the  wavefunction  for  a diatomic  molecule  can  be 
written  as 


^ vjmA'  X.  <#>)  (29) 

where  (r,  R)  is  an  electronic  wavefunction  for  state  i at  fixed  inter- 
nuclear  separation  R,  (R)  is  a vibrational  wavefunction  for  level  v and 

vJnA  (9,  x5  0)  refers  to  the  rotational  state  specified  by  electronic  angular 
momentum  A,  total  angular  momentum  J and  magnetic  quantum  number  M.  The  rep- 
resentation is  in  a coordinate  system  related  to  a space-fixed  system  by  the 
Eulerian  angles  (0,  0) . The  transition  moment  Mjj_  can  be  written,  using 

the  wavefunction  given  by  Eq.  (29),  as 
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The  subscripts  e,  v and  r refer  to  the  electronic,  vibrational  and  rotational 
wavefunct ions  and  and  M11  are  the  electronic  and  nuclear  electric  dipole 
moments,  respectively.  Integration  over  the  electronic  wavefunct ion,  in  the 
Born-Oppenheimer  approximation,  causes  the  contribution  of  the  nuclear  moment 
Mn  to  vanish  for  i ^ j.  The  electronic  dipole  moment  can  be  written  (Refs.  51 
and  52)  in  the  form 


Me  = -Ier  (3D 

where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordinates 
rk  refer  to  a molecule-fixed  system  and  $ (0,  x,  0)  is  a group  rotation  ten- 
sor whose  elements  are  the  direction  cosines  related  to  the  Eulerian  rotation 
angles  (9,  x,  0).  Using  bracket  notation,  Eqs . (30)  and  (31)  can  be  combined 
to  yield  for  the  transition  moment 
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(32) 


The  matrix  elements  (J*  A'  M'  I*  (e,  X,  0)|j"  A"  M")  determine  the  group  selec- 
tion rules  for  an  allowed  transition  and  have  been  evaluated  for  many  types  of 
transitions  (Refs.  53,  54  and  55).  Summing  Eq.  (32)  over  the  degenerate  mag- 
netic quantum  numbers  M'  and  M"  we  have  from  Eq.  (28) 
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where  is  the  Honl-London  factor  (Refs.  56  and  57)  and 
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is  the  band  strength  for  the  transition.  Combining  Eqs.  (25),  (27)  and  (33), 
we  have  for  the  intensity  of  a single  emitting  line  from  upper  level  n: 


Inm  = I mV"j"  3 NJ' 
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[^Emv"j*J  smv"j"A" 


ri4c3u>n 
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where  Nj  is  the  number  density  in  the  upper  rotational  state  J and  o>n  = 
(2-6{j  A)  (2S*  + l)  is  the  electronic  degeneracy.  Taking  an  average  value  of 
Em  v"  J”  *>or  t*ie  whole  band,  Eq.  (35)  can  be  summed  to  yield  the  total  inten- 
sity in  the  (v’,  v")  band: 
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where  Nv’  = j , j,  is  the  total  number  density  in  the  upper  vibrational  level 
v’  and  where  we  make  use  of  the  group  summation  property 


Comparing  Eq.  (27)  and  (36),  we  have  for  the  Einstein  spontaneous  transition 
coefficient  of  the  band  (v’,  v") 
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Similarly,  the  lifetime  of  an  upper  vibrational  level  v’  of  state  n can  be 
written 
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where  the  summat ion  runs  over  all  v"  for  each  lower  state  m.  Equation  (3 c) 
can  be  cast  in  the  computational  form 
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wh  -re  £e£  '!/'  and  p_“  y”  are  in  atomic  units.  It  is  also  often  convenient  to 
relate  the  transition  probability  to  the  number  of  dispersion  electrons  d~d 
to  explain  the  emission  strength  classically.  This  number,  the  f-numb . • or 
oscillator  strength  for  emission,  is  given  by 
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The  inverse  process  of  absorption  is  related  to  th*“  above  development 
through  the  Einstein  5 coefficient.  Corresponding  to  Eq.  (27)»  we  hav- 
single  line  in  absorption 


T - T nV'J' 

I mn  “ ^mv"j"  ' 


/ 


K (1/ 


di/  - AEmnNmBmn 


line  (v'Vj"j') 


(42) 


where  K(u)  is  the  absorption  coefficient  of  a beam  of  photons  of  frequency  u 
and 
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(43) 


is  the  Einstein  absorption  coefficient  for  a single  line.  Summing  over  all 
lines  in  the  band  (v",  v’),  assuming  an  average  band  frequency,  we  obtain 
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where  Nv"  = j"  j " is  the  total  number  density  in  the  lower  vibrational  state 
v".  Corresponding  to  Eqs.  (40)  and  (4l)  we  can  define  an  f-number  or  oscil- 
lator strength  for  absorption  as 
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In  computational  form,  Eq.  (45)  becomes 
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where  &E^  and  p”  are  in  atomic  units.  Combining  Eqs.  (38)  and  (4l)  and 
comparing  with  Eq.  (45),  we  see  that  the  absorption  and  emission  f-numbers  are 
related  by 


f mn,v"v' 


(47) 


Some  caution  must  be  observed  in  the  use  of  f-numbers  given  either  by  Eq.  (4l) 
or  (45)  since  both  band  f-numbers  and  system  f-numbers  are  defined  in  the  lit- 
erature. The  confusion  arises  from  the  several  possible  band  averaging  schemes 
that  can  be  identified. 

An  integrated  absorption  coefficient  (density  corrected)  can  be  defined 
from  Eq.  (44)  as 
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where  the  exponential  factor  corrects  for  stimulated  emission.  Equation  (48) 
can  be  written  in  terms  of  the  absorption  f-number  as 
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Using  h c/k  = 1.43880  cm-K° , we  obtain  a computational  formula  for  the 
integrated  absorption  coefficient  as  Sv»5  v<  (cm-2  . atm-1)  = 


2 3795  x I07 
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The  total  integrated  absorption  is  found  from 


stotal=X  X SV«V' 

V"  V'  (51) 

where,  under  normal  temperature  conditions,  only  the  first  few  fundamentals 
and  overtones  contribute  to  the  summations. 

The  developments  given  above  are  rigorous  for  band  systems  where  an 
average  band  frequency  can  be  meaningfully  defined.  Further  approximations, 
however,  are  often  made.  For  example,  the  electronic  component  of  the  dipole 
transition  moment  can  be  defined  as 


<||-Ierji>  (52) 

This  quantity  is  often  a slowly  varying  function  of  R and  an  average  value  can 
sometimes  be  chosen.  Equation  (34)  can  then  be  written  approximately  in  fac- 
tored form  as 


Pmv"  - Qv'v'Z  kji  <R)|2 


(53) 


where  qv»  v»,  the  square  of  the  vibrational  overlap  integral,  is  called  the 
Franck -Condon  factor.  ^jl  ls  evaluated  at  some  mean  value  of  the  inter- 

nuclear  separation  R.  In  addition,  it  is  sometimes  possible  to  account  for  a 
weak  R-dependence  in  by  a Taylor  series  expansion  of  this  quantity  about 
some  reference  value,  /P^q,  usually  referred  to  the  (0,  0)  band.  We  have 
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Substituting  into  Eq.  (3*0  and  integrating  yields 
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where 
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is  the  R-centroid  for  the  transition  and 
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is  the  R -centroid.  Note  that  this  last  term  differs  (to  second  order)  from 
the  square  of  the  R-centroid.  An  alternate  procedure  can  be  developed  by 
evaluating  Eq.  (52)  at  each  R-centroid,  Ryi  «.  Then 


Pmv"  **  qv/v"£|/?j<(RvV,)| 
1 I J 


(58) 


Equation  (58)  assumes  that  the  vibrational  wavefunction  product  w , t.  behaves 
like  a delta  function  upon  integration. 


Vv  vv; = ^(r-  Rv'v")  <^y' I v'/> 


(59) 


The  range  of  validity  of  Eq.  (58)  is  therefore  questionable,  particularly  for 
band  systems  with  bad  overlap  conditions  such  as  oxygen  Schumann-Runge . The 
range  of  validity  of  the  R-centroid  approximation  has  been  examined  by  Frazer 
(Ref.  58). 
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The  final  step  in  calculating  transition  probabilities  is  the  determination 
of  #■;?/,  *he  electronic  dipole  transition  moment,  for  the  entire  range  of 
inter:,  iclear  separations,  R,  reached  in  the  vibrational  levels  to  be  consid- 
ered. This  can  be  expressed  in  terms  of  the  expansion  of  Eq.  (4)  as 


*,,(«)=!  c'*cj,  <^(R)|  Me|*„(R)> 

ii  1/  ' 


where  cu^  and  c*  are  coefficients  for  and  ;e]_ , respectively. 

An  analysis  similar  to  that  yielding  Eqs . (10)  and  (ll)  gives 
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The  spatial  integral  in  Eq.  (6l)  reduces  to  one-electron  integrals  equivalent 
to  overlap  integrals,  and  the  evaluation  of  Eq.  (6l)  can  be  carried  out  by  the 
same  computer  programs  used  for  Eq.  (ll).  Programs  for  evaluating  /?.^(R)  in 
Eq.  (60)  have  been  developed  at  UTRC  and  examples  of  their  application  have 
appeared  in  the  literature  (Ref.  59) • 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have 
a strong  R-dependence  and  R-centroid  or  other  approximations  will  be  invalid. 

A direct  evaluation  of  Eq.  (34)  would  therefore  be  required  using  the  fully- 
coupled  system  of  electronic  and  vibrational  wavefunctions  to  properly  account 
for  the  source  of  the  band  perturbations. 
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SECTION  IV 


DISCUSSION  OF  RESULTS 


theoru*  leal  research  conducted  under  this  program  was  concerned  with 
1-  * r : . i ■ tructure  and  radiative  transition  probabilities  for  10,  NO?  and 
his  research  was  carried  out  using  quantum  mechanical  methods  and  corn- 
put*  r codes  which  were  developed  at  this  Center.  Both  ab  initio  methods, 

: pl./irv  c nfi.yuration  interaction  (Cl),  and  density  functional  methods  (X^), 

■ Llized.  The  results  obtained  for  the  molecular  systems  under  study  are 
summarized  separately  below.  . 

MO 


An  extensive  analysis  of  the  dipole  moment  function  for  the  ground  state 
of  VO  was  carried  out.  We  have  examined  in  detail  the  recent  work  of  Goldman 
and  Schmidt  (Ref.  60),  Chandraiah  and  Cho  (Ref.  6l),  and  the  theoretical 
studies  of  Billingsley  (Refs.  62  and  63).  In  addition  we  have  reviewed  our 
valence-bond  calculations  carried  out  for  both  the  separated  atom  range  of 
interaction  of  U+0  and  for  the  united  atom  limit  of  the  state  of  phosphorus. 
Computer  least  square  fits  to  all  the  available  experimental  and  theoretical 
data  relating  to  the  dipole  moment  were  carried  out  and  our  resultant  best  fit 
is  shown  in  Fig.  1.  This  curve  is  qualitatively  similar  to  that  reported  pre- 
viously (Ref.  64)  but  differs  significantly  in  several  respects.  First,  the 
permanent  dipole  moment  of  TO  is  negative  corresponding  to  N"0+  polarity.  This 
is  .now  well  established  both  from  theoretical  studies  and  comparisons  with  the 
experimental  data  for  CO.  Second,  the  short  range  portion  of  p(R)  shows  a 
minimum  between  .5  - -6  X,  well  up  on  the  repulsive  part  of  the  potential  curve 
but  not  yet  at  united  atom  behavior  for  the  system.  Third,  the  long  range 
behavior  of  p (R ) has  now  been  calculated  with  high  accuracy  optimized  -wavefunc- 
tions . This  synthesized  curve  for  the  dipole  moment  is  the  best  that  can  be 
obtained  using  all  of  the  data  available  at  the  present  time. 

All  data,  including  experimental  integrated  absorption  coefficients,  were 
finally  reduced  to  fit  a dipole  moment  function  of  the  form: 

^ ( r ) L a (r-re)J  (62) 

J=0  J 


A computer  reduction  was  carried  out  using  a least -squares  fit  to  all  data 
weighted  by  the  given  error  analysis  (or  estimated  errors  in  the  case  of  the 
calculated  dipole  moment  functions).  A minimum  error  residual  was  found  with 
the  following  function: 

N0(r)  = -0.166  + 2.24  ( r-re ) - 1.342  (r-re)2  - 1.703  (r-rj3  + 1.097  (r-r ^ (63) 
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where  re  is  the  equilibrium  internuclear  separation  in  angstroms  and  p,(r)  is 
in  Debyes. 

This  function  can  be  compared  with  the  "experimental"  dipole  moment 
function  for  CO  found  by  Toth,  et  al  (Ref.  65): 


, CO 
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(64) 


The  functions  are  clearly  similar  both  with  regard  to  the  sign  of  the  coeffi- 
cients and  their  relative  magnitudes. 

Using  the  dipole  moment  function  represented  by  Eq.  (63)  and  accurate  RKR 
vibrational-rotational  wavefunctions , the  integrated  band  absorption  coeffi- 
cients have  been  calculated  for  the  fundamental,  first,  and  second  overtones 
of  NO  as  a function  of  temperature.  For  a transition  from  lower  state  1 to 
upper  state  u in  absorption  we  have 
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where  f^u  is  the  absorption  f -number.  In  computational  form  we  have 
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The  total,  absorption  is  given  by 


S(T)  =£Sj|u(T) 
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where  the  first  term  represents  the  total  absorption  of  the  fundamental  band, 
the  second  term  represents  the  first  overtone  absorption,  etc.  The  calculated 
absorption  is  given  in  Table  1 which  indicates  a near  constancy  of  the  inten- 
sity of  the  fundamental  band  up  to  5 000°K . This  is  in  agreement  with  the 
recent  experimental  data  of  Konkov  and  Vorontsov  (Ref.  66).  Table  2 lists 
fundamental  oscillator  strengths  and  Table  3 lists  the  calculated  Einstein 
coefficients  for  the  lowest  10  vibrational  levels. 
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NOg 

SCF  calculations  were  performed  for  the  ground  X ^A^state  of  NO2  at  the 
experimental  equilibrium  bond  length  of  1.1934  A and  at  bond  angles  of  134.1° 
(exp.)  and  180°.  The  results  are  given  in  Table  4 where  we  list  the  occupied 
molecular  orbitals  at  these  two  bond  angles.  The  bent  molecule  has  a lower 
total  energy  in  the  ground  state  and  the  lowest  occupied  orbitals  correspond 
to  those  predicted  from  ab  initio  calculations.  The  density  functional  method 
thus  successfully  predicts  the  bent  character  of  the  ground  state  and  the  cor- 
rect orbital  occupancy.  Table  4 also  lists  the  energies  of  the  low-lying 
4b2,  7aj_  and  2b^  states. 

Additional  calculations  were  carried  out  for  the  first  excited  A state 
of  NC^.  The  calculations  were  performed  at  the  experimental  ground  state  bond 
length  and  at  bond  angles  of  134.1°  and  180°.  The  results  are  given  in  Table  5 
which  lists  the  occupied  molecular  orbitals.  The  linear  configuration  is 
favored  in  this  excited  state. 

From  this  data  we  calculate  a vertical  excitation  energy  for  X 2A]_  -•  A ^B]_ 
of  2.43  eV.  This  is  in  good  agreement  with  experimental  estimates  and  the  best 
ab  initio  calculation  for  this  system  (Ref.  67).  More  detailed  calculations  of 
the  NO2  hypersurfaces,  using  the  density  functional  method,  are  clearly  indi- 
cated as  a result  of  the  success  of  these  initial  studies. 

TiO 

SCF  calculations  were  performed  for  the  d state  of  TiO  at  four 
internuclear  separations  (R  = 2.5,  3.0>  3*5  and  4.0  bohrs),  using  an  optimized 
SCF  orbital  basis  obtained  at  R = 3.0.  This  orbital  basis  is  shown  in  Table 
6.  The  MO's  obtained  from  these  calculations  were  then  used  as  input  trans- 
formation vectors  to  convert  the  original  atomic  orbital  basis  to  MO  form. 

The  lowest  eleven  (11)  MO's  (la,  2a,  3a,  1tt+,  lrr~,  4a,  5a,  2n  + , 2tt- , 6a,  7a) 
were  taken  as  doubly  occupied  and  a Cl  was  performed  over  the  rest  of  the  MO's 
for  the  X^a,  a^A  and  d-^E*  states  of  TiO.  This  resulted  in  a wavefunction 
expansion  which  consisted  of  135  configurations  for  X^a , 102  configurations  for 
a^A  and  88  configurations  for  d^E+.  A full  Cl  was  not  performed  due  to  the 
number  of  configurations  involved.  Those  configurations  selected  were  con- 
structed from  the  lowest  9 sigma  MD's,  the  lowest  4 pi  MO's  and  a single  delta 
MO.  This  corresponds  to  the  energetically  more  compacted  orbitals. 

Expectation  values  for  the  dipole  moment  p(R)  were  calculated  for  the 
X^a , a^  and  d1E+ states.  The  results  are  shown  in  Fig.  2. 

The  X^a,  a-^-A  and  d^E+  vibrational  states  of  TiO  can  be  characterized  by 
the  spectroscopic  data  (Ref.  68)  tabulated  in  Table  7.  A RKR  analysis  of  this 
data  yields  the  potential  curves  shown  in  Fig.  3.  These  potentials  were  then 
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used  to  generate  wavefunctions  for  the  lowest  twenty  (20)  vibrational  levels 
of  the  x3& , a^A  and  d^2+states  of  TiO.  The  calculated  AG(v)  and  experimental 
AG(v)  are  compared  in  Tables  8,  9,  and  10. 

Using  these  vibrational  wavefunctions  in  conjunction  with  the  dipole 
moment  function  shown  in  Fig.  2,  the  vibrational  oscillator  strengths  were 
calculated.  These  are  shown  in  Tables  11,  12  and  13  for  several  of  the  low- 
lying  fundamental  and  first  and  second  overtones.  The  population  densities 
as  a function  of  temperature  are  shown  in  Table  lU  for  those  vibrational  levels 
which  were  significantly  populated.  The  three  states  (X^a,  a^A  and  d^2  4 ) 
were  weighted,  assuming  a Boltzmann  population  distribution,  for  calculation 
of  integrated  absorption  coefficients.  These  absorption  data  are  shown  in 
Table  15. 

Density  functional  studies  were  undertaken  on  the  calculation  of  the 
electronic  structure  of  TiO  and  TiO+  using  the  SCF-Xq,  scattered  wave  method. 

An  efficient  version  of  the  Xq,  program  has  been  developed  in  this  Center  for 
carrying  out  these  studies.  In  addition,  a separate  computer  program  has  been 
developed  for  implementing  the  scattered-wave  approach  for  atoms  and  atomic 
ions  incorporating  both  a nonrelativistic  option  and  a relativistic  option 
using  a modified  Eriet  Hamiltonian.  The  results  of  our  SCF-Xq,  calculations 
for  the  constituent  atoms  and  ions  are  given  in  Table  16.  For  these  calcula- 
tions, the  exchange  parameter  a was  optimized  to  satisfy  the  virial  theorem. 
This  procedure  yields  a total  energy  which  is  in  very  good  agreement  with  the 
exact  Hartree  Fock  results.  Examination  of  the  energies  given  in  Table  l6 
shows  that  s-d  hybridization  will  be  very  important  in  TiO  owing  to  the  near 
degeneracy  of  the  orbital  energies  in  Ti.  We  also  see  that  the  Ti+0“structure 
should  be  dominant  for  low-lying  electronic  states  of  TiO  which  are  compati- 
ble with  this  symmetry.  These  observations  are  borne  out  by  the  molecular 
SCF-Xq.  results  for  TiO  (R  = 3-0  a.u.)  given  in  Table  17.  The  dominant  atomic 
component  of  each  molecular  orbital  is  given  in  Table  17.  For  both  (16^) 
and(l6  5°)  configurations,  the  molecule  is  predicted  to  be  highly  ionic.  The 
closeness  of  the  MO  eigenvalues  for  these  different  configurations  indicates 
that  there  should  be  low-lying  2 and  A states  for  this  molecule.  The  experi- 
mental data  shown  in  Table  17  are  consistent  with  the  predicted  SCF-Xq,  results. 
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TABLE  1 


Total  Integrated  Absorption  Coefficients 
For  Nitric  Oxide  (X  ^H) 


Absorption 

Coefficient, 

p _■] 

cm  atm 

Temperature 

(°K) 

Fundamental 

First 

Overtone 

Second 

Overtone 

Third 

Overtone 

Total 

100. 

122.725 

2.101 

.029 

.000 

124.855 

273.16 

122.725 

2.102 

.029 

.000 

124.856 

300. 

122.724 

2.102 

.029 

.000 

124.855 

500. 

122.689 

2.120 

.03  0 

.000 

124.839 

1000. 

122.149 

2.402 

.036 

.000 

124.587 

1500. 

121.131 

2.932 

.051 

.001 

124.115 

2000. 

119.889 

3.575 

.073 

.001 

123.538 

2500. 

118.530 

4.274 

.102 

.002 

122.908 

3000. 

117.098 

5.006 

.140 

.004 

122.248 

4000. 

114.087 

6.519 

.242 

.009 

120.857 

5000. 

110.922 

8.040 

.376 

.018 

119.356 

8000. 

100.644 

11.927 

.848 

.055 

113.474 

10000. 

93.513 

13.552 

1.108 

0 

00 

108.251 
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TABLE  2 


Calculated  Oscillator  Strengths  (fv,vn) 
for  the  X ^tt  - X ^ir  System  of  NO 


f ** 

' V 

0 

1 

2 

3 

4 

5 

0 

1 

5.158-06 

2 

8.631-03 

9.982-06 

3 

1.224-09 

2.627-07 

1.447-05 

4 

1.224-11 

4.912-09 

5.222-07 

1.363-05 

5 

5.076-14 

7.748-11 

1.241-03 

8.660-07 

2.245-05 

6 

1.084-14 

7.943-14 

2.676-10 

2.500-08 

1.293-06 

2.593-05 

7 

7.117-14 

7.134-15 

1.310-12 

6.700-10 

4.422-08 

1.799-06 

3 

9.898-14 

9.419-14 

1.203-13 

8.250-12 

1.376-09 

7.182-08 

9 

1.267-14 

8.275-14 

1.117-13 

4.086-13 

2.952-11 

2.526-09 

10 

1.630-14 

1.135-14 

1.384-13 

8.905-15 

2.499-13 

6.758-11 

■'/v" 

6 

7 

3 

9 

10 

0 

1 

2 

3 

4 


9 

10 


2 . 909-05 

2.382-06 

1.093-07 

4.369-09 


3.191-05 

3.040-06 

1.581-07 


3.441-05 

3.769-06 


3.658-05 
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TABLE  3 


! 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Calculated  Einstein  Coefficients  (A 

p p V v"  ' 

the  X ^rr  - X it  System  of  NO  ’ 


for 


v'/v” 

0 

1 

2 

3 

4 

5 

0 

1 

1.210  + 01 

2 

8.165  - 01 

2.273  +01 

3 

2.509  - 02 

2.357  +00 

3.196+01 

4 

4.464  - 04 

9.762  - 02 

4.542  +00 

3.988  +01 

5 

2.302  - 06 

2.696  - 03 

2.392  - 01 

7.300  +00 

4.656  +01 

6 

8.482  - 07 

4.251  - 06 

9.022  - 03 

4.667-  01 

1.056  +01 

5.209  + 01 

7 

7.465  - 06 

5.450  - 07 

6.792  - 05 

2.188  - 02 

7.994  - 01 

1.422  +01 

3 

1.335  - 05 

9.573  - 06 

8.841  - 06 

4.143  - o4 

4.352  - 02 

1.257+00 

9 

2.129  - 06 

1.031  - 05 

1.099  - 05 

2.907  - 05 

1.435  - 03 

7.730-  02 

10 

3.430-  06 

1.027  - 06 

2.382  - 05 

8.483  - 07 

1.720-  05 

3.177-03 

v’/v” 

6 

7 

8 

9 

10 

5.655  +01 
1.822  +01 
1.850  +00 
1.293  - 01 


6.002  +01 
2.249  +01 
2.588  + 00 


6.256  + 01 
2.694  +01 


6.426  +01 


4o 


TAB  Li;  4 


Total  and 

/.  o rbital  Energies 

for  the  Oround  X A State 

of  NO 

1 

2 

rt  i a 

: * it  an  ry 

rbital  Energy 

(P.yd.  ) 

134.1 

180° 

la  ( 1 c ) 

p 

-37.9329 

-37.9957 

b2  (lau} 

d 

-37.9329 

-37.9957 

2a,  (2a  ) 
X g 

2 

-28.8059 

-28.8329 

3a  (3a  ) 
1 g 

2 

- 2.1963 

- 2.2209 

2b  (2a  ) 
a 

2 

- 1.9703 

- 2.0685 

4a  ( 4a  ) 
1 3 

2 

- 1.2915 

- 1.3266 

5a.,  (Itt  ) 
1 u 

2 

- 1.0217 

- 1.0690 

3b2  (3au) 

2 

- 1.1126 

- 1.1571 

lb,  (in  ) 
1 u 

2 

- 1.0396 

- 1.0690 

4b0  (in  ) 
g 

2 

- O.7898 

- 0.8902 

la  ( Itt  ) 
2 g 

2 

- 0.8392 

- 0.8902 

6a,  (2tt  ) 
1 u 

1 

- 0.6177 

- 0.4741 

2b.  (2n  ) 

x U 

0 

- 0.4410 

- 0.4741 

7a-,  (5a  ) 

g 

0 

- 0.08492 

- 

5b2  (4au) 

0 

- 

- 

Total  Energy (Rydbergs ) -403. 1783 

= 2.2552  a. u.  (1.1934  A) 

n-o 


4l 


TABLE  5 


~ 2 

Total  and  Orbital  Energies  for  the  First  Excited  A B]_  State  of  NO2 


Orbital 

Occupancy 

Orbital  Energy 

(Ryd. ) 

13l .lu 

180° 

lal  (^g) 

2 

-37.9180 

-37.9957 

lb2  (lau) 

2 

-37.9180 

-37.9957 

2al  (2CTg) 

2 

-28.7221 

-28.8329 

Sa-L  ( 3ag) 

2 

- 2.1755 

- 2.2209 

2b2  (2au) 

2 

- 1.9539 

- 2.0685 

la.  (la  ) 
1 g 

2 

- 1.2720 

- 1.3266 

5al  (lTTu} 

2 

- 1.0012 

- 1.0690 

3b 2 (3cu) 

2 

- 1.091+5 

- 1.1571 

lb,  (In  ) 
1 u 

2 

- 1.0192 

- 1.0690 

l+b?  (in  ) 
d g 

2 

- 0.7781 

- 0.8902 

la9  (lrr  ) 

2 g 

2 

- 0.8281 

- 0.8902 

6al  <2ttu) 

0 

- 0.5968 

- O.I7I1 

2bx  (2nu) 

1 

- 0.1ll3 

- O.I7I1 

(5cg) 

0 

- 0.07899 

- 

5b  (la  ) 
2 u 

0 

- 

- 

Total  Energy  (Rydbergs)  -102.9990 

R A = 2.2552  a.  u.  (1.1931+  A) 

N-0 


h2 


TABLE  6 


MOLECULAR  ORBITAL  BASIS  FOR  TiO 


Molecular  Orbital 


Screening  Parameters  for  Atomic  Orbitals 


1 c 

Is  (Ti) 

2cr 

2s  (Ti) 

3 a 

Is  (0) 

It r 

2Pt^»-  (Ti) 

1 *c 

2pa  (Ti) 

5a 

3s  (Ti) 

2-rr  »_ 

3p7r+’“  (Ti) 

6a 

3pO"  (Ti) 

7a 

2s  (0) 

3/- 

2p7r+>“  (0) 

8a 

2pa  (0) 

9 a 

l*s  (Ti) 

1*7^*- 

3d*''’"  (Ti) 

10  a 

3da  (Ti) 

11a 

l*pa  (Ti ) 

i8+ 

3dS+  (Ti) 

21. 1*1*09 
7.6883 
7.6579 
9.0321* 
9.0321* 
3.6777 
3.3679 
3.3679 
2.2187 
2.0885 
2.2211* 
1.3973 
2.2035 
2.8021 
1.0852 
2.7000 


ABLE  8 


VALUES  OF  AG(v)  FOR  THE  X 3 A STATE  OF  TiO 
v AG(v) , cm  3‘ 


Calculated 

0 

999.08 

1 

989.86 

2 

980.61 

3 

971.12 

1 

962.20 

5 

952.99 

6 

913.77 

7 

931.55 

8 

925.33 

9 

916.11 

10 

906.89 

11 

897.67 

12 

888.16 

13 

879-21 

ll 

870.03 

15 

860.81 

16 

851.65 

17 

812.18 

18 

833.29 

19 

821.00 

20 

Experimental 
999.18 
989.96 
980. 7l 
971.52 
962.30 
953.08 

9U3.86 
931.61 
925.12 
916.20 
906.98 
897.76 
888. 5k 
879.32 
870.10 


TABLE  10 

Values  of AG(v)  for  the  d 1£+  State  of  TiO 


V 


AG(v),  cm-1 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


Calc ilated 

1014 . 50 
1005 . 30 

996 . 10 

986.90 

977.70 

968.50 

959.30 
950.1C 

940.91 

931.71 

922 . 51 

913.31 

904.11 

894.92 
885.73 

876.55 

867.39 

858.22 

849.03 

839-71 

47 


Experimental 

1014.60 

1005-40 

996.20 

987.00 
977-80 

968.60 

959-40 

950.20 

941.00 
931-80 

922.60 

913-40 

904.20 

895.00 
885.80 

876.60 
867.40 

858.20 

849.00 

4 39-00 


UmA 


uMftfll 


TABLE  11 


CALCULATED 

OSCILLATOR  STRENGTHS 

FOR  THE  d 1£+  STATE 

OF  TiO 

Oscillator  Strengths 

Vibrational  Level 

Fundamental 

First  Overtone 

Second  Overtone 

Vi" 

V Vl«" 

V vi+2" 

fVi"  vi+3" 

0 

5.5755-06 

1.5557-07 

1.2650-08 

1 

1.1168-05 

3.531+9-07 

1.8265-08 

2 

1.6366-05 

5. 9581-07 

3.1+928-08 

3 

2.11+55-05 

8.979I-07 

5.1+611-08 

1+ 

2.6323-05 

1.2172-06 

7.7966-08 

5 

3.1000-05 

1.5781-06 

1.1186-07 

6 

3.5571-05 

1.9583-06 

1.5172-07 

7 

1+.  0003-05 

2.3302-06 

2.081+1-07 

8 

1.1*2*28-05 

2.6893-06 

2.8511-07 

9 

1+.  8901-05 

3.0266-06 

3.6192-07 

10 

5.3289-05 

3.1+039-06 

1*.  1975-07 

11 

5.71*1*1-05 

3.8751-06 

1+ . 5915-07 

12 

6.1301+-05 

1*.  1+077-06 

1+.  8956-07 

13 

6.1+902-05 

5.0183-06 

5.0251-07 

lit 

6.8135-05 

5.7662-06 

1*.  9720-07 

15 

7.0919-05 

6.61+52-06 

1+.  7758-07 

16 

7.3209-05 

7.6988-06 

1+.1+1+81+-07 

17 

7.1*933-05 

8.9351+-06 

1+.  2136-07 

18 

7.6193-05 

1.0302-05 

— 

19 

7.6999-05 

— 

— 

i 


1+8 


T’A^LE  12 


CALCULATED  OSCILLATOR  STRENGTHS  FOR  THE  a 1A  STATE  OF  TiO 


Vibrational  Level 
Vi" 

Oscillator  Strengths 

Fundamental  First  Overtone  Second  Overtone 

V n.i'  V vit2"  f-'i" 

0 

1.1207-05 

2.3931-08 

2.1780-08 

1 

2.0939-05 

9.0683-08 

3.61+72-08 

2 

3-0069-05 

2.6922-07 

1+.  581+2-08 

3 

3.8578-05 

5.6097-07 

6.1+891-08 

1+ 

1+. 6191-05 

9.893U-07 

7.9385-08 

5 

5.308U-05 

1.6057-06 

9.0295-08 

6 

5.9200-05 

2.3986-06 

1.0172-07 

7 

6. 1+U59-05 

3.3893-06 

1.0573-07 

8 

6.8966-05 

1+. 621+3-06 

1.0831-07 

9 

7.2611-05 

6.1313 -06 

1.1630-07 

10 

7.5251-05 

7.9721+-06 

1.2327-07 

11 

7.6932-05 

1.0190-05 

1.2821-07 

12 

7.76LL-05 

1.2819-05 

1.3755-07 

13 

7.7273-05 

1.5962-05 

1.5093-07 

lU 

7.5802-05 

1.9733-05 

1.6950-07 

15 

7.3198-05 

2.1+213-05 

1.8893-07 

16 

6.9525-05 

2.9338-05 

1.8572-07 

17 

6.5103-05 

3.1+958-05 

1.5075-07 

18 

6.0153-05 

1*.  1027-05 

— 

19 

5.1+689-05 

— 

— 

>+9 


'"ABLE  13 


CALCULATED  OSCILLATOR  STRENGTHS  FOR  THE  X STATE  OF  TiO 


Vibrational  Level 
vi" 

Oscillator  Strength 
Fundamental  First  Overtone 

f*  II  II  f , , 

vi  vi+1  Vi"  vi+2 

s 

Second  Overtone 

f„  » v 11 

vi  vi+3 

0 

6.2728-06 

2. 3090-08 

3.2107-09 

1 

1.2008-05 

7.5950-08 

1*.  9075-09 

2 

1. 71*78-05 

1.8398-07 

5.20148-09 

3 

2.2688-05 

3.U581-07 

6J4569-O9 

1* 

2.75V-05 

5.6595-07 

6.1(1473-09 

5 

3.211U-05 

8.538SW07 

l*.  6627-09 

6 

3.6II0-05 

I.I8I48-O5 

1.7328-09 

T 

^.0595-05 

1.5353-06 

8.1605-13 

8 

14.1*660-05 

1.9303-06 

8. 71*33-10 

9 

14.8)409-05 

2.1*578-06 

2.8862-10 

10 

5.1I48I4-05 

3.19I4I4-O6 

3.0370-10 

11 

5.3870-05 

l*.  1391-06 

3.5653-09 

12 

5.5589-05 

5.361I4-06 

1.6738-08 

13 

5.6373-05 

6.9961-06 

1*.  8901-08 

lU 

5.61142-05 

9.1306-06 

I.106I4-07 

15 

5 - 1+827-05 

1.1917-05 

2.0735-07 

l6 

5.21437-05 

1.5308-05 

2.9961+-07 

17 

l*.  9356-05 

1.9176-05 

3.6017-07 

18 

I4.57U6-O5 

2.352U-05 

— 

19 

I4.16I4I4-05 

— 

— 
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.'ABLE  14 


f 


POPULATION  DENSITIES  FOR  THE  X 3 A , a XA  AND  d 1£  + STATES  OF  TiO 


Vibrational 


Population  Density 


T (°K) 

Level 

x 3A 

a 2a 

d *2* 

273.16 

0 

9.79513-01 

1.53049-02 

9.10073-07 

1 

5.07401-03 

7.92817-05 

4.34702-09 

300. 

0 

9.71734-01 

1.99653-02 

2.66473-06 

1 

8.06081-03 

1.65618-04 

2.05310-08 

2 

6.98899-05 

1.43596-06 

1.65337-10 

3 

6.33366-07 

1. 30131-08 

1.39166-12 

500. 

0 

8.87709-01 

5.55980-02 

1.99469-04 

1 

5.00691-02 

3.13587-03 

1.07629-05 

2 

2.89996-03 

1.81627-04 

5.96353-07 

3 

1.72479-04 

1.08025-05 

3.39314-08 

4 

1.05342-05 

6.59766-07 

1.98254-09 

1000. 

0 

6.61852-01 

9.56300-02 

4.05030-03 

1 

1.57185-01 

2.27114-02 

9. 40835-04 

2 

3.78287-02 

5.46581-03 

2.21463-04 

3 

9.22557-03 

1.33299-03 

5.28263-05 

4 

2.27996-03 

3.29427-04 

1.27691-05 

5 

5.70980-04 

8.25000-05 

3.12775-06 

6 

1.44903-04 

2.09368-05 

7-76364-07 

7 

3.72643-05 

5.38427-06 

1.95281-07 

8 

9.71117-06 

1.40315-06 

4.97754-08 

1500. 

0 

5.08061-01 

9.69984-02 

9.35441-03 

1 

1.94842-01 

3.71990-02 

3.53482-03 

2 

7.53858-02 

1.43926-02 

1.34759-03 

3 

2.94264-02 

5.61805-03 

5.18313-04 

4 

1.15885-02 

2.21246-03 

2.01125-04 

5 

4.60422-03 

8.79031-04 

7.87372-05 

6 

1.84555-03 

3.52350-04 

3.10982-05 

7 

7.46342-04 

1.42491-04 

1.23917-05 

8 

3.04502-04 

5.81351-05 

4.98159-06 

9 

1.25338-04 

2.39294-05 

2.02044-06 

10 

5.20495-05 

9.937''  I— 06 

8.26731-07 

11 

2.18068-05 

4.16;-.",  '-06 

3.41290-07 
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TABLE  l4  (Cont'd) 


T (°K) 

Vibrational 

Level 

X 3 A 

Population  Density 
a 1A 

d xz  + 

2000. 

0 

4.07242-01 

8.93735-02 

1.30059-02 

1 

1. 98462-01 

4.35545-02 

6.26835-03 

2 

9-73605-02 

2.13668-02 

3.04122-03 

3 

4.80805-02 

1.05518-02 

1.48533-03 

4 

2.39021-02 

5.24555-03 

7. 30259-04 

5 

1.19614-02 

2.62506-03 

3.61420-04 

6 

6.02574-03 

1.32241-03 

1.80065-04 

7 

3.05576-03 

6.70618-04 

9.03080-05 

8 

1.55994-03 

3.42345-04 

4.55936-05 

9 

8.01637-04 

1.75927-04 

2.31719-05 

10 

4.14694-04 

9.10089-05 

1.18550-05 

11 

2.15953-04 

4.73931-05 

6.10549-06 

12 

1.13206-04 

2.48443-05 

3.16534-06 

13 

5.97396-05 

1.31105-05 

1.65197-06 

lU 

3.17348-05 

6.96452-06 

8.67886-07 

15 

1.69703-05 

3.72430-06 

4.58991-07 

2500. 

0 

3.  37858-01 

8.06111-02 

1.50153-02 

1 

1.90105-01 

4.53581-02 

8.37423-03 

2 

1.07537-01 

2.56578-02 

4.69526-03 

3 

6.11543-02 

.1.45911-02 

2.64655-03 

4 

3.49623-02 

8.34183-03 

1.49970-03 

5 

2.00945-02 

4.79445-03 

8.54343-04 

6 

1.16108-02 

2.77026-03 

4.89289-04 

7 

6.74446-03 

1.60919-03 

2.81710-04 

8 

3.93857-03 

9-39721-04 

1.6 3059-04 

9 

2.31225-03 

5.51690-04 

9.48837-05 

10 

1.36469-03 

3.25608-04 

5.55064-  : 

11 

8.0/728-04 

1.93197-04 

3.2'  436- 

12 

4.83002-04 

1.15242-04 

1.93000-05 

13 

2.89642-04 

6.91071-05 

1.14716-05 

14 

1.74615-04 

4.16621-05 

15 

1.05829-04 

2.52501-05 

4.11780-  • 

16 

6.44808-05 

1.53848-05 

!.  48662-  • 

17 

3.94068-05 

0.42373-06 

1.5  #62-04 

3000. 

0 

2.87715-01 

7.25815-02 

1.59381- 

1 

1.78174-01 

4.49479-02 

9.70743-03 

2 

1.10828-01 

0.70534-02 

i . )4  > V - 3 • 

3 

6.92425-02 

1 . 74677-02 

3.75167-03 
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Population  Density 


TABLE  14  (Cont'd) 


Vibrational 

Level 

x 3 A 

a XA 

d XI  + 

4 

4.34528-02 

1.09618-02 

2.33701-03 

5 

2.73894-02 

6.90949-03 

1.46224-03 

6 

1.73407-02 

4.37453-03 

9.18959-04 

7 

1.10274-02 

2.78187-03 

5.80088-04 

8 

7.04367-03 

1.77690-03 

3.67801-04 

9 

4.51903-03 

1.14001-03 

2.34235-04 

10 

2.91214-03 

7.34642-04 

1.49834-04 

11 

1.88495-03 

4.75514-04 

9.62699-05 

12 

1.22548-03 

3.09151-04 

6.21284-05 

13 

8.00268-04 

2.01883-04 

4.02727-05 

l4 

5.24909-04 

1.32418-04 

2.62212-05 

15 

3.45822-04 

8.72402-05 

1.71480-05 

l6 

2.28845-04 

5.77306-05 

1.12641-05 

17 

1.52108-04 

3.83721-05 

7.43185-05 

18 

1.01551-04 

2.56180-05 

4.92515-06 

19 

6.80978-05 

1.71789-05 

3.27840-06 

TABLE  15 


INTEGRATED  ABSORPTION  COEFFICIENTS 

FOR  TiO 

T (°K) 

STotal  (cm'2  atrr_1) 

STotal  (cm~‘ 

273.16 

145.86 

1U5.86 

300. 

146.77 

133. 64 

500. 

155.  42 

84.91 

1000. 

161.25 

44.05 

1500. 

160. 23v 

18.25 

2000. 

157.39 

21.49 

2500. 

154.69 

16.90 

3000. 

153.51 

13.98 

S _ £ j_+]_  where  the  summation  extends  over  all  populated 

Total  i fundamental  bands . 

S*Total  = STotal  x = density  corrected  absorption 

L 'J  coefficient. 


+The  absorption  coefficient  should  normally  increase  as  temperature  in- 
creases. The  decrease  shown  from  1500°K  to  3000°K  is  probably  due  to 
small  errors  in  the  relative  dipole  moment  functions  for  X , 
and  d d£+  states  of  TiO. 
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TABLE  16:  SCF-X  RESULTS  FOR  ATOMIC  SPECIES 
a 


SPECIES 

ORBITAL 

OCCUPATION 

ORBITAL  ENERGY 

2 2 

Ti  (4d  6s  - non  rel. ) 

Is 

2s 

2p 

3s 

2 

0 

-355.3260 

- 39.0396 

- 32.7064 

- 4.4908 

Total  Energy: 
-I696.6824  Ryd. 

c. 

6 

2 

3p 

6 

- 2.8168 

3d 

2 

- O.2968 

4s 

2 

- 0.2880 

Ti  (4d^  6s  - rel.) 

Is 

2 

-356.6870 

Total  Energy: 

2s* 

2 

- 39.2258 

-1699.9234  Ryd. 

2p 

2 

- 32.7804 

2p 

4 

- 32.3576 

3s* 

2 

- 4.4802 

3p 

2 

- 2.7890 

3P* 

4 

- 2.7402 

3d 

2 

- 0.2446 

3d 

0 

- 0.2414 

4s 

2 

- O.2732 

+ 2 

Ti  (4d  6s  - non  rel.) 

Is 

2 

-355.8690 

Total  Energy: 

2s 

2 

- 39-5728 

-I696.I808  Ryd. 

2p 

6 

- 33-2406 

3s 

2 

- 5.0270 

3p 

6 

- 3.3526 

3d 

2 

- 0.8198 

4s 

1 

- 0.7244 

4 

0 (2p  - non  rel. ) 

Is 

2 

- 37.8192 

Total  Energy: 

2s 

2 

- 1.7304 

-149.5210  Ryd. 

2p 

4 

- 0.6556 

(Ryd. 
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TABLE  17:  SCF-X  ORBITAL  ENERGIES  BOR  TiO 

' a 


ORBITAL 

OCCUPATION 

ORBITAL  ENERGY  (Ryd.) 

2 

TiO  (16  ) 

Ti  (is) 

2 

-355-3052 

1’^Z+ , ]T  (ave) 

Ti  (2s) 

2 

- 39-0229 

0 (is) 

2 

- 38.1172 

R=3.0  a.u. 

Ti  (2p) 

6 

- 32.6894 

la(  Ti-3s) 

2 

- 4.5036 

2a( Ti-3p) 

2 

- 2.8475 

3a(0-2s) 

2 

- 1.1362 

4a(0-2p) 

2 

- 0.4239 

5a(Ti-4s) 

0 

- 

ln(Ti-3p) 

4 

- 2.8357 

2n(0-2p) 

4 

- 0.5093 

3n(Ti-3d) 

0 

- 

16(Ti-3d) 

2 

- 0.3617 

TiO  (165a) 

Ti  (is) 

2 

-355.2238 

T ( ave . ) 

Ti  (2s) 

2 

- 38.9438 

0 (Is) 

2 

- 38.3109 

R=3.0  a.u. 

Ti  (2p) 

6 

- 32.6097 

la( Ti-3s) 

2 

- 4.4467 

2ct( Ti-3p) 

2 

- 2.7952 

3ct(0-2s) 

2 

- 1.1849 

4ct(0  -2p) 

2 

- 0.4101 

5a(Ti-4s) 

1 

- O.316I 

1tt(T1-3p) 

4 

- 2.7818 

2n(0  -2p) 

4 

- 0.5196 

3n(Ti-3<i) 

0 

- 

16(Ti-3d) 

1 

- 0.3259 
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ABSTRACT 


Calculations  have  been  performed  for  the  kinetics  of  the  dissociative- 
recombination  of  an  electron  and  the  positive  ion  of  nitric  oxide.  Elec- 
tronic wavefunctions  have  been  constructed  for  selected  dissociating  states 
of  NO  and  expectation  values  of  the  electronic  energy  and  electronic  transi- 
tion matrix  elements  have  been  calculated.  The  electronic  wavefunctions  for 
NO  were  analyzed  to  determine  those  states  responsible  for  dissociative- 

o -f  p 2 

recombination.  States  of  , ^tt  and  symmetry  are  found  to  be  the  most 
important  channels  for  this  reaction.  Calculations  are  presented  for  the 
dxssociation-recombination  cross-sections  as  a function  of  electron  energy 
and  vibrational  state  of  the  molecular  ion.  These  data  are  also  reduced 
to  kinetic  rate  constants  as  a function  of  the  electron  temperature  and  ion 
vibrational  temperature,  assuming  independent  Boltzmann  distribution  func- 
tions. We  find  a calculated  temperature  dependence  of  T^-'  J for  the  rate 
coefficient  with  branching  primarily  into  N + 0 ^P. 
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J.  S.  Wagner,  Physics  Dept 

Govt  Fiscal  Relations  & Patent  Off 
275  Admin  Building  AG-70 
University  of  Washington 
Seattle  WA  98195 
Attn:  R.  Geballe 

K.  C.  Clark 
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! ■ ; ' iit-ractorc 

(Cont ' J 

Harvard  'nivrrsity 

its  pheric  Science: 

Fieri?-  nail 
Card  rl  1 • MA  . 13- 
Att . M Etroy 


. Air  a Circle 
(edi>ri,  MA  01730 
Attn:  D.  Hansen 

M.  P.  Shuler 
H.  Stewart 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 
Attn:  H.  Wolfhard 

E.  Bauer 
J.  Ben/s ton 
P.  A.  Selwyn 

IHT  Corporation 

P.  0.  Box  1087 
San  Diego,  CA  92111 
Attn:  J.  A.  Rutherford 

R.  H.  Neynaber 
D.  A.  Vroom 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel  MD  20810 

Attn:  Document  Librarian 


Kaman  Sciences  Corporation 
P.0.  Box  7463 

Colorado  Springs,  CO  80933 
Attn:  B.  J.  Bittner 

P.  Jesson 
F.  H.  Shelton 

LIT  Research  Institute 
10  West  35th  Street 
Chicago,  111.  60616 

Attn:  Technical  Library 


Lockheed  Missiles  and  Space  Company 
•151  Hanover  Street 
Palo  Alto,  CA  94304 

Attn:  M.  Walt  Dept  52-10 

J.  B.  Regan  d/52-12 
B.  M.  McCormac  Dept  52 -l4 
T.  James 

J.  L.  Kulander  Dept  52-14 
J.  Kumer 

R.  D.  Sears  Dept  52- l4 
R.  G.  Johnson  Dept  52-12 
R.  C.  Gunton  Dept  52-14 


Maryland,  University  of 
Administration  Building 
College  Park,  MD  21030 
Attn:  T.  J.  Rosenberg 

R.  J.  Munn  Dept  of  Chem 
J.  Vanderslice  Chem  Dept 
Dept  of  Physics  H.  Griem 

Maxwell  Laboratories,  Inc. 

'■244  Balboa  Avenue 
San  Diego,  CA  92123 
Attn:  V.  Fargo 

McDonnell  Dougles  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Attr : J.  Moule 

A.  D.  Goedeke 
W.  Olsen 


Photometries,  Inc. 

442  Marrett  Road 
Lexington,  MA  02173 
Attn:  I.  L.  Kofsky 

Physical  Dynamics  Inc. 
P.0.  Box  1069 
Berkeley,  CA  94701 
Attn:  A.  Thompson 

J.  B.  Workman 


Department  of  Defense  Contractors 
(Cont'd) 

Physical  Sciences,  Inc. 

30  Commerce  Way 
Woburn,  MA  01801 


RScD  Associates 

1815  N.  Ft.  Myer  Drive 

] 1th  Floor 

Arlington,  VA  22209 
Attn:  H.  J.  Mitchell 

J.  M.  Rosengren 


Attn:  R.  L.  Taylor 

K.  Wray 

Physics  International  Company 

2700  Merced  Street 

San  Leandro,  CA  94577 

Attn:  DOC  CON  FOR  Tech  Lib 


Pittsburgh,  University  of 
Of  the  Comwlth  Sys  of  Higher  Educ 
Cathedral  of  Learning 
Pittsburgh,  PA  15213 
Attn:  M.  A.  Biondi 

E.  Gerjudy 

F.  Kaufman 
W.  L.  Fite 

Princeton  Univ.  The  Trustees  of 
Forres tal  Campus  Library 
Box  36 

Princeton  University 
Princeton,  NJ  08540 
Attn:  A.  J.  Kelly 


R&D  Associates 
P.0.  Box  9695 

Marina  Del  Re.y,  CA  90291 
Attn:  B.  Gabbard 

R.  Latter 
F.  Gilmore 

C.  H.  Humphrey 

D.  Dee 

R.  P.  Turco 
H.  A.  Dry 
R.  G.  Lindgren 
R.  E.  Lelevier 
A.  L.  Latter 


Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  90406 
Attn:  C.  Crain 

Rice  University 
Department  of  Space  Science 
Houston,  TX  77058 

Attn:  R.  F.  Stebbings 

Midwest  Research  Institute 
425  Volker  Boulevard 
Kansas  City,  MO  64110 

Attn:  Physics  Sec  T.  H.  Milne 
Physics  Sec  F.  T.  Greene 

Minnesota,  University  of 
2030  University  Avenue,  S.E. 
Minneapolis,  MN  554l4 
(No  Classified  to  this  Address) 
Attn:  J.  R.  Winkler 

Minnesota,  University  of 
Morris  Campus 
Morris , MN  56267 
Attn:  M.  Hersch 


Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  CA  93101 
Attn:  R.  Hendrick 

C.  L.  Longmire 
R.  Kilb 


D.  Archer 

M.  Scheibe 
R.  N.  Wilson 
D . Sowle 

D.  Sappenfield 
W.  F.  Crevier 

N.  J.  Carron 
R.  Bogus ch 
P.  Fischer 
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Department  of  Defense  Contractors 
(Cont 'd) 

Mitre  Corporation,  The 
P.0.  Box  208 
Bedford,  MA  01730 
Attn:  S.  A.  Morin  M/S 

J.  N.  Freedman 

National  Academy  of  Sciences 
Attn:  Nat'l  Materials  Advisory  Board 
2101  Constitution  Avenue 
Washington,  D.C.  204l8 
Attn:  W.  C.  Bartley 

J.  R.  Sievers 

E.  R.  Dyer  (2  copies) 

New  York  at  Buffalo,  State  Univ.  of 
Buffalo,  NY  14214 
Attn:  G.  0 Brink 

Pacific-Sierra  Research  Corp. 

1456  Cloverfield  Blvd 
Santa  Monica,  CA  90404 
Attn:  E.  C.  Field  Jr 


Riverside  Research  Institute 
80  West  End  Avenue 
New  York,  NY  10023 
Attn:  G.  Glaser 

J.  B.  Minkoff 
D.  H.  Koppel 

Science  Applications,  Inc. 

P.0.  Box  2351 
La  Jolla,  CA  92038 
Attn:  D.  A.  Hamlin 

B.  Myers 
R.  W.  Lowen 

Science  Applications,  Inc. 
Huntsville  Division 
2109  W.  Clinton  Avenue 
Suite  700 

Huntsville,  AL  35805 
Attn:  N.  R.  Byrn 

Professor  Chalmers  F.  Sechrist 
155  Electrical  Engineering  Bldg. 
University  of  Illinois 
Urbana,  IL  6l801 


Pennsylvania  State  University 
Industrial  Security  Office 
Room  55  Old  Main  Bldg 
University  Park,  PA  16802 
Attn:  L.  Hale 

J.  S.  Nisbet 


Rice  University,  William  Marsh 
P.0.  Box  2692 
Houston,  TX  77001 
.Attn:  Industrial  Sec  Super 

J.  Chamberlain 


Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Attn:  W.  G.  Chestnut  (2  copies) 

B.  R.  Gasten 
A.  L.  Whitson 
A.  M.  Peterson 
R.  White 

R.  L.  Leadabrand 
J.  T.  Moseley 
T.  G.  Slanger 
D.  J.  Hildenbrand 
L.  L.  Cobb 
J.  R.  Peterson 
J.  Lomax 
R.  D.  Hake  Jr. 

F.  T.  Smith  (UNCL  Only) 

G.  Black 
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Department  of  Defense  Contractors 
(Cont'd) 

Stanford  Research  Institute 
306  Wynn  Drive,  N.W. 

Huntsville,  AL  35805 
Attn:  M.  Morgan 

Systems,  Science  and  Software,  Inc. 
P.0.  Box  1620 
La  Jolla,  CA  92038 
Attn:  R.  C.  Sklarew 

H.  B.  Levine 

Technology  International  Corporation 
75  Wiggins  Avenue 
Bedford,  MA  01730 
Attn:  W.  P.  Boquist 

TRW  Systems  Group 

One  Space  Park 

Redondo  Beach,  CA  90278 

Attn:  Tech  Infor  Center/S-1930 

R.  Watson  Rl/1096 
J.  F.  Friichtenicht  Rl-1196 


United  Technologies  Research  Center 
East  Hartford,  CT  06108 

Attn:  B,  H.  Bullis  & H.  H.  Michels 

Utah  State  University 
Logan,  UT  84321 

Attn:  W.  M.  More  Chem  Dept 

D.  Burt 

C.  Wyatt 

D.  Baker 
K.  Baker 


Visidyne,  Inc. 

19  Third  Avenue 
North  West  Industrial  Rirk 
Burlington,  MA  01803 
Attn:  R.  Waltz 

J.  W.  Carpenter 
T.  C.  Degges 
0.  Manley 


Wayne  State  University 
1064  Mackenzie  Hall 
Detroit,  MI  48202 
Attn:  P.  K.  Rol 

R.  H.  Kummler 

Wayne  State  University 
Dept,  of  Physics 
Detroit,  MI  48202 
Attn:  W.  E.  Kauppila 

Yale  University 
New  Haven,  CT  06520 

Attn:  Engrg  Dept.  G.  J.  Schultz 


Denver,  University  of 
Colorado  Seminary 
Denver  Research  Institute 
P.  0.  Box  10127 
Denver,  CO  80210 

Attn:  SEC  OFF  FOR  D.  Muccray 
SEC  OFF  FOR  V.  Zyl 

Howard  University 
Department  of  Chemistry 
Washington,  D.  C.  20059 
Attn:  W.  Jackson 

Johns  Hopkins  University 
34th  & Charles  Streets 
Baltimore,  I©  21218 
Attn:  J.  J.  Kauftaan 

Southern  California,  Univ.  of 
University  Park 
Los  Angeles,  CA  90007 
Attn:  S.  W.  Benson 

Stanford  Research  Institute 
l6ll  North  Kent  Street 
Arlington,  VA  22209 
Attn:  W.  W.  Berning 
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Continued  from  Department  of  the  Army 


Commander 

Harry  Diamond  Laboratories 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 
Attn:  M.  H.  Weiner 

: DRYD-NP  (2  copies) 


Commander 

Trasana 

White  Sands  Missile  Range  NM  88002 
Attn:  R.  E.  Dekinder,  Jr. 

Director 

U.  S.  Army  Ballistic  Res.  Labs. 

Aberdeen  Proving  Ground,  MD  21005 
Attn:  M.  D.  Kregel 

DRXBR-AM  J.  S.  Batten 
J.  C.  Mester 
Tech  Lib  E.  Baicy 
J.  A.  Vanderhoff 
DRXBR-AM  D.  Snider 
DRXBR-CA  T.  J.  O'Brien 
J.  Heimerl  AMXRD-BSP 
I.  L.  Chidsey 
DRXBR-AM  G.  E.  Keller 


Chief  C-E  Services  Division 
U.  S.  Army  Communications  CMD 
Pentagon  Rm  13269 
Washington,  D.  C.  20310 
Attn:  CEAO 

Commander 

US  Army  Electronics  Command 
Fort  Monmouth,  N.J.  07703 

Attn:  DRSEL-PL-ENV  H.  A.  Bomke 


